Вектор Лапласа — Рунге — Ленца

В этой статье векторы выделены жирным шрифтом, а их абсолютные величины — курсивом, например, $|\mathbf {A} |=A$.

Ве́ктор Лапла́са — Ру́нге — Ле́нца (вектор Лапласа, вектор Рунге — Ленца и вектор Ленца) — вектор, который используется для описания формы и ориентации орбиты, по которой одно небесное тело обращается вокруг другого (например, орбиты, по которой планета обращается вокруг звезды). В случае с двумя телами, взаимодействие которых описывается законом всемирного тяготения Ньютона, вектор Лапласа — Рунге — Ленца представляет собой интеграл движения, то есть его направление и величина постоянны независимо от точки орбиты, в которой они вычисляются; говорят, что вектор Лапласа — Рунге — Ленца сохраняется при гравитационном взаимодействии двух тел. Это утверждение можно обобщить на любую задачу с двумя телами, взаимодействующими посредством центральной силы, которая изменяется обратно пропорционально квадрату расстояния между ними. Такая задача называется Кеплеровой задачей.

Например, такой потенциал возникает при рассмотрении классических орбит (без учёта квантования) в задаче о движении отрицательно заряженного электрона в электрическом поле положительно заряженного ядра. Если вектор Лапласа — Рунге — Ленца задан, то форма относительного движения тел может быть получена из простых геометрических соображений, с использованием законов сохранения этого вектора и энергии.

Согласно принципу соответствия, у вектора Лапласа — Рунге — Ленца имеется квантовый аналог, который был использован в первом выводе спектра атома водорода ещё до открытия уравнения Шрёдингера.

В задаче Кеплера существует особенность: конец вектора импульса p всегда движется по окружности. Из-за расположения этих кругов, для заданной полной энергии E, задача Кеплера математически эквивалентна задаче о частице, свободно перемещающейся в четырёхмерной сфере $S_{3}$. Согласно этой математической аналогии, сохраняющийся вектор Лапласа — Рунге — Ленца соответствует дополнительным компонентам углового момента в четырёхмерном пространстве.

Вектор Лапласа — Рунге — Ленца также известен как вектор Лапласа, вектор Рунге — Ленца и вектор Ленца, хотя ни один из этих учёных не был его первооткрывателем. Вектор Лапласа — Рунге — Ленца переоткрывался несколько раз. Он также эквивалентен безразмерному вектору эксцентриситета в небесной механике. Для него также нет никакого общепринятого обозначения, хотя обычно используется символ A. Для различных обобщений вектора Лапласа — Рунге — Ленца, которые будут определены ниже, используется символ ${\mathcal {A}}$.

Контекст

Для одиночной частицы, движущейся под воздействием любой консервативной центральной силы, существуют по крайней мере четыре интеграла движения (сохраняющиеся величины): полная энергия E и три компоненты вектора углового момента $\mathbf {L}$. Орбита частицы лежит в плоскости, определяемой начальным импульсом частицы p или скоростью v и её радиус-вектором r (рис. 1). Эта плоскость перпендикулярна постоянному вектору $\mathbf {L}$, что можно выразить математически с помощью скалярного произведения $\mathbf {r} \cdot \mathbf {L} =0$.

Как указано ниже, вектор Лапласа — Рунге — Ленца A всегда находится в плоскости движения, то есть равенство $\mathbf {A} \cdot \mathbf {L} =0$ выполняется для любой центральной силы. Он также является постоянным только для силы, зависящей обратно пропорционально квадрату расстояния. Если центральная сила приближённо зависит от обратного квадрата расстояния, вектор $\mathbf {A}$ является почти постоянным по длине, но медленно вращается. Для большинства центральных сил этот вектор $\mathbf {A}$ не постоянен и изменяет как длину, так и направление. Обобщённый сохраняющийся вектор Лапласа — Рунге — Ленца ${\mathcal {A}}$ может быть определён для всех центральных сил, но этот вектор представляет собой сложную функцию положения и обычно не выражается аналитически в элементарных или специальных функциях.

История

Вектор Лапласа — Рунге — Ленца A является сохраняющейся величиной в задаче Кеплера и полезен при описании астрономических орбит, например движения планеты вокруг Солнца. Однако он никогда не был широко известен среди физиков, возможно, потому что он менее интуитивно понятен, чем импульс и угловой момент. Вектор Лапласа — Рунге — Ленца независимо открывали несколько раз за прошедшие три столетия. Якоб Герман был первым, кто показал, что вектор A сохраняется для специального случая центральной силы, зависящей обратно пропорционально квадрату расстояния, и нашёл его связь с эксцентриситетом эллиптической орбиты. Работа Германа была обобщена до её современной формы Иоганном Бернулли в 1710 году. В свою очередь, Пьер-Симон Лаплас в конце XVIII столетия заново открыл сохранение вектора $\mathbf {A}$, доказав это аналитически, а не геометрически, как его предшественники.

В середине XIX века Уильям Гамильтон получил эквивалент вектора эксцентриситета, определённый ниже, и использовал его, чтобы показать, что конец вектора импульса $\mathbf {p}$ движется по кругу под действием центральной силы, зависящей обратно пропорционально квадрату расстояния (рис. 3). В начале XX столетия Уиллард Гиббс нашёл тот же самый вектор с помощью векторного анализа. Вывод Гиббса использовал Карл Рунге в популярном немецком учебнике по векторам в качестве примера, на который ссылался Вильгельм Ленц в своей статье о квантовомеханическом рассмотрении атома водорода.

В 1926 году этот вектор применил Вольфганг Паули для вывода спектра атома водорода, используя современную матричную квантовую механику, а не уравнение Шрёдингера. После публикации Паули вектор стал известен как вектор Рунге — Ленца.

Математическое определение

Рис. 1. Вектор Лапласа — Рунге — Ленца A (показанный красным цветом) в четырёх точках (обозначенных 1, 2, 3 и 4) на эллиптической орбите связанной точечной частицы, движущейся под действием центральной силы, которая зависит обратно пропорционально квадрату расстояния. Маленький чёрный круг обозначает центр притяжения. От него начинаются радиус-векторы (выделены чёрным цветом), направленные в точки 1, 2, 3 и 4. Вектор углового момента L направлен перпендикулярно орбите. Компланарные векторы $\mathbf {p} \times \mathbf {L}$, $(mk/r)\mathbf {r}$ и A изображены синим, зелёным и красным цветами, соответственно; эти переменные определены ниже. Вектор A является постоянным по направлению и величине

Для одиночной частицы, движущейся под воздействием центральной силы, зависящей обратно пропорционально квадрату расстояния и описываемой уравнением $\mathbf {F} (\mathbf {r} )={\frac {-k}{r^{2}}}\mathbf {\hat {r}}$, вектор Лапласа — Рунге — Ленца A определён математически формулой

$\mathbf {A} =\mathbf {p} \times \mathbf {L} -mk\mathbf {\hat {r}} ,$

где $m$ — масса точечной частицы, движущейся под воздействием центральной силы, $\mathbf {p}$ — вектор импульса, $\mathbf {L} =\mathbf {r} \times \mathbf {p}$ — вектор углового момента, $k$ — параметр, описывающий величину центральной силы, $\mathbf {\hat {r}}$ — единичный вектор, то есть $\mathbf {\hat {r}} ={\frac {\mathbf {r} }{r}}$, где $\mathbf {r}$ — радиус-вектор положения частицы, и $r$ — его длина.

Поскольку предполагается, что сила консервативная, то полная энергия системы $E$ сохраняется

$E={\frac {p^{2}}{2m}}-{\frac {k}{r}}={\frac {1}{2}}mv^{2}-{\frac {k}{r}}.$

Из центральности силы следует, что вектор углового момента L также сохраняется и определяет плоскость, в которой частица движется. Вектор Лапласа — Рунге — Ленца A перпендикулярен вектору углового момента L и, таким образом, находится в плоскости орбиты. Уравнение $\mathbf {A} \cdot \mathbf {L} =0$ верно, потому что векторы $\mathbf {p} \times \mathbf {L}$ и r перпендикулярны L.

Это определение вектора Лапласа — Рунге — Ленца A применимо для одиночной точечной частицы с массой $m$, движущейся в стационарном (не зависящем от времени) потенциале. Кроме того, это определение может быть применено к задаче двух тел, такой как задача Кеплера, если заменить $m$ на приведённую массу этих двух тел и r на вектор между ними.

Круговой годограф импульса

Рис. 2. Конец вектора импульса p (показанный синим цветом) двигается по кругу, когда частица совершает движение по эллипсу. Четыре помеченные точки соответствуют точкам на рис. 1. Центр круга находится на оси $y$ в точке $A/L$ (показан пурпурным), с радиусом $mk/L$ (показан зелёным). Угол $\eta$ определяет эксцентриситет $e$ эллиптической орбиты ($\cos \eta =e$). Из теоремы о вписанном угле для круга следует, что $\eta$ является также углом между любой точкой на окружности и двумя точками пересечения окружности с осью $p_{x}$, $p_{x}=\pm p_{0}$

Сохранение вектора Лапласа — Рунге — Ленца A и вектора углового момента L используется в доказательстве того, что конец вектора импульса $\mathbf {p}$ движется по окружности под действием центральной силы, зависящей обратно пропорционально квадрату расстояния. Вычисляя векторное произведение A и L, получается уравнение для p

$L^{2}\mathbf {p} =\mathbf {L} \times \mathbf {A} -mk{\hat {\mathbf {r} }}\times \mathbf {L} .$

Направляя вектор $\mathbf {L}$ вдоль оси $z$, а главную полуось — вдоль оси $x$, получаем уравнение

$p_{x}^{2}+(p_{y}-A/L)^{2}=(mk/L)^{2}\,.$

Другими словами, конец вектора импульса p движется по окружности радиуса $mk/L$, центр которой расположен в точке с координатами $(0,A/L)$. Эксцентриситет $e$ равен косинусу угла $\eta$, показанного на рис. 2. Для краткости вводится переменная $p_{0}={\sqrt {2m|E|}}$. Круговой годограф полезен для описания симметрии задачи Кеплера.

Интегралы движения и суперинтегрируемость

Семь скалярных величин — энергия $E$ и компоненты векторов Лапласа — Рунге — Ленца A и момента импульса L — связаны двумя соотношениями. Для векторов выполняется условие ортогональности $\mathbf {A} \cdot \mathbf {L} =0$, а энергия входит в выражение для квадрата длины вектора Лапласа — Рунге — Ленца, полученного выше $A^{2}=m^{2}k^{2}+2mEL^{2}$. Тогда существует пять независимых сохраняющихся величин, или интегралов движения. Это совместимо с шестью начальными условиями (начальное положение частицы и её скорость являются векторами с тремя компонентами), которые определяют орбиту частицы, так как начальное время не определено интегралами движения. Поскольку величину A (и эксцентриситет $e$ орбиты) можно определить из полного углового момента $L$ и энергии E, то утверждается, что только направление $\mathbf {A}$ сохраняется независимо. Кроме того, вектор A должен быть перпендикулярным L — это приводит к одной дополнительной сохраняющейся величине.

Механическая система с d степенями свободы может обладать максимум $2d-1$ интегралами движения, поскольку имеется $2d$ начальных условий, а начальное время не может быть определено из интегралов движения. Система с более чем $d$ интегралами движения называется суперинтегрируемой, а система с $2d-1$ интегралами называется максимально суперинтегрируемой. Поскольку решение уравнения Гамильтона — Якоби в одной системе координат может привести только к d интегралам движения, то переменные должны разделяться для суперинтегрируемых систем в больше чем одной системе координат. Проблема Кеплера максимально суперинтегрируема, так как она имеет три степени свободы ($d=3$) и пять независимых интегралов движения; переменные в уравнении Гамильтона — Якоби разделяются в сферических координатах и параболических координатах, как описано ниже. Максимально суперинтегрируемые системы могут быть квантованы с использованием только коммутационных соотношений, как показано ниже.

Уравнение Гамильтона — Якоби в параболических координатах

Постоянство вектора Лапласа — Рунге — Ленца можно вывести, используя уравнение Гамильтона — Якоби в параболических координатах $(\xi ,\eta )$, которые определяются следующим образом

$\xi =r+x\,,$

$\eta =r-x\,,$

где $r={\sqrt {x^{2}+y^{2}}}$ — радиус в плоскости орбиты. Обратное преобразование этих координат запишется в виде

$x={\frac {1}{2}}(\xi -\eta )\,,$

$y={\sqrt {\xi \eta }}\,.$

Разделение переменных в уравнении Гамильтона — Якоби в этих координатах даёт два эквивалентных уравнения

$2\xi p_{\xi }^{2}-mk-mE\xi =-\beta \,,$

$2\eta p_{\eta }^{2}-mk-mE\eta =\beta \,,$

где $\beta$ — интеграл движения. Посредством вычитания этих уравнений и выражения в терминах декартовых координат импульса $p_{x}$ и $p_{y}$ можно показать, что $\beta$ эквивалентен вектору Лапласа — Рунге — Ленца

$\beta =p_{y}(xp_{y}-yp_{x})-mk{\frac {x}{r}}=A_{x}\,.$

Этот подход Гамильтона — Якоби может использоваться, чтобы вывести сохраняющийся обобщённый вектор Лапласа — Рунге — Ленца ${\mathcal {A}}$ в присутствии электрического поля E

${\mathcal {A}}=\mathbf {A} +{\frac {mq}{2}}\left[(\mathbf {r} \times \mathbf {E} )\times \mathbf {r} \right]\,,$

где $q$ — заряд обращающейся частицы.

Альтернативная формулировка

В отличие от импульса p и углового момента L, у вектора Лапласа — Рунге — Ленца нет общепринятого определения. В научной литературе используются несколько различных множителей и символов. Самое общее определение даётся выше, но другое определение возникает после деления на постоянную $mk$, чтобы получить безразмерный сохраняющийся вектор эксцентриситета

$\mathbf {e} ={\frac {1}{mk}}(\mathbf {p} \times \mathbf {L} )-\mathbf {\hat {r}} ={\frac {m}{k}}(\mathbf {v} \times \mathbf {r} \times \mathbf {v} )-\mathbf {\hat {r}} \,,$

где v — вектор скорости. Направление этого масштабированного вектора $\mathbf {e}$ совпадает с направлением A, и его амплитуда равна эксцентриситету орбиты. Мы получим другие определения, если поделить $\mathbf {A}$ на $m$:

$\mathbf {M} =\mathbf {v} \times \mathbf {L} -k\mathbf {\hat {r}}$

или на $p_{0}$:

$\mathbf {D} ={\frac {\mathbf {A} }{p_{0}}}={\frac {1}{\sqrt {2m|E|}}}\{\mathbf {p} \times \mathbf {L} -mk\mathbf {\hat {r}} \}\,,$

который имеет ту же размерность, что и угловой момент (вектор L). В редких случаях знак вектора Лапласа — Рунге — Ленца может быть изменён на противоположный. Другие общие символы для вектора Лапласа — Рунге — Ленца включают a, R, F, J и V. Однако выбор множителя и символа для вектора Лапласа — Рунге — Ленца не влияет на его сохранение.

Рис. 3. Вектор углового момента L, вектор Лапласа — Рунге — Ленца A и вектор Гамильтона, бинормаль B, являются взаимно перпендикулярными; A и B указывают соответственно на большую и на малую полуоси эллиптической орбиты в задаче Кеплера

Альтернативный сохраняющийся вектор: бинормаль — вектор B изучен Уильямом Гамильтоном

$\mathbf {B} =\mathbf {p} -\left({\frac {mk}{L^{2}r}}\right)(\mathbf {L} \times \mathbf {r} ),$

который сохраняется и указывает вдоль малой полуоси эллипса. Вектор Лапласа — Рунге — Ленца $\mathbf {A} =\mathbf {B} \times \mathbf {L}$ является векторным произведением B и L (рис. 3). Вектор B обозначен как бинормаль, так как он перпендикулярен как A, так и L. Подобно вектору Лапласа — Рунге — Ленца, вектор бинормали можно определить с различными множителями.

Два сохраняющиеся вектора A и B можно объединить в сохраняющийся двухэлементный тензор

$\mathbf {W} =\alpha \mathbf {A} \otimes \mathbf {A} +\beta \mathbf {B} \otimes \mathbf {B} \,,$

где $\otimes$ обозначает тензорное произведение, а $\alpha$ и $\beta$ — произвольные множители. Записанное в компонентной записи, это уравнение читается так

$W_{ij}=\alpha A_{i}A_{j}+\beta B_{i}B_{j}\,.$

Векторы $\mathbf {A}$ и $\mathbf {B}$ ортогональны друг другу, и их можно представить как главные оси сохраняющегося тензора $\mathbf {W}$, то есть как его собственные вектора. $\mathbf {W}$ перпендикулярен $\mathbf {L}$:

$\mathbf {L} \cdot \mathbf {W} =\alpha (\mathbf {L} \cdot \mathbf {A} )\mathbf {A} +\beta (\mathbf {L} \cdot \mathbf {B} )\mathbf {B} =0\,,$

поскольку $\mathbf {A}$ и $\mathbf {B}$ перпендикулярны, то $\mathbf {L} \cdot \mathbf {A} =\mathbf {L} \cdot \mathbf {B} =0$.

Вывод орбит Кеплера

Рис. 4. Упрощённая версия рис. 1. Определяется угол $\theta$ между $\scriptstyle \mathbf {A}$ и $\scriptstyle \mathbf {r}$ в одной точке орбиты.

Форму и ориентацию орбиты в задаче Кеплера, зная вектор Лапласа — Рунге — Ленца $\mathbf {A}$, можно определить следующим образом. Рассмотрим скалярное произведение векторов $\mathbf {A}$ и $\mathbf {r}$ (положения планеты):

$\mathbf {A} \cdot \mathbf {r} =Ar\cos \theta =\mathbf {r} \cdot (\mathbf {p} \times \mathbf {L} )-mkr,$

где $\theta$ является углом между $\mathbf {r}$ и $\mathbf {A}$ (рис. 4). Поменяем порядок множителей в смешанном произведении $\mathbf {r} \cdot (\mathbf {p} \times \mathbf {L} )=\mathbf {L} \cdot (\mathbf {r} \times \mathbf {p} )=\mathbf {L} \cdot \mathbf {L} =L^{2}$, и при помощи несложных преобразований получим определение для конического сечения:

${\frac {1}{r}}={\frac {mk}{L^{2}}}\left(1+{\frac {A}{mk}}\cos \theta \right)$

с эксцентриситетом $e$, заданным по формуле:

$e={\frac {A}{mk}}={\frac {|\mathbf {A} |}{mk}}.$

Приходим к выражению квадрата модуля вектора A в виде

$A^{2}=m^{2}k^{2}+2mEL^{2}\,,$

которое можно переписать, используя эксцентриситет орбиты

$e^{2}-1={\frac {2L^{2}}{mk^{2}}}E\,.$

Таким образом, если энергия отрицательна, что соответствует связанным орбитам, эксцентриситет меньше, чем единица, и орбита имеет форму эллипса. Наоборот, если энергия положительна (несвязанные орбиты, также называемые орбитами рассеяния), эксцентриситет больше, чем единица, и орбита — гипербола. Наконец, если энергия точно равна нулю, эксцентриситет — единица, и орбита — парабола. Во всех случаях вектор $\mathbf {A}$ направлен вдоль оси симметрии конического сечения и указывает на точку самого близкого положения точечной частицы от начала координат (перицентр).

Сохранение под действием силы, обратно пропорциональной квадрату расстояния

Сила $\mathbf {F}$, действующая на частицу, предполагается центральной. Поэтому

$\mathbf {F} ={\frac {d\mathbf {p} }{dt}}=f(r){\frac {\mathbf {r} }{r}}=f(r)\mathbf {\hat {r}}$

для некоторой функции $f(r)$ радиуса $r$. Поскольку угловой момент $\mathbf {L} =\mathbf {r} \times \mathbf {p}$ сохраняется под действием центральных сил, то ${\frac {d}{dt}}\mathbf {L} =0$ и

${\frac {d}{dt}}(\mathbf {p} \times \mathbf {L} )={\frac {d\mathbf {p} }{dt}}\times \mathbf {L} =f(r)\mathbf {\hat {r}} \times \left(\mathbf {r} \times m{\frac {d\mathbf {r} }{dt}}\right)=f(r){\frac {m}{r}}\left[\mathbf {r} \left(\mathbf {r} \cdot {\frac {d\mathbf {r} }{dt}}\right)-r^{2}{\frac {d\mathbf {r} }{dt}}\right],$

где импульс записан в виде $\mathbf {p} =m{\frac {d\mathbf {r} }{dt}}$, и двойное векторное произведение упростилось с помощью формулы Лагранжа

$\mathbf {r} \times \left(\mathbf {r} \times {\frac {d\mathbf {r} }{dt}}\right)=\mathbf {r} \left(\mathbf {r} \cdot {\frac {d\mathbf {r} }{dt}}\right)-r^{2}{\frac {d\mathbf {r} }{dt}}.$

Тождество

${\frac {d}{dt}}(\mathbf {r} \cdot \mathbf {r} )=2\mathbf {r} \cdot {\frac {d\mathbf {r} }{dt}}={\frac {d}{dt}}(r^{2})=2r{\frac {dr}{dt}}$

приводит к уравнению

${\frac {d}{dt}}(\mathbf {p} \times \mathbf {L} )=-mf(r)r^{2}\left[{\frac {1}{r}}{\frac {d\mathbf {r} }{dt}}-{\frac {\mathbf {r} }{r^{2}}}{\frac {dr}{dt}}\right]=-mf(r)r^{2}{\frac {d}{dt}}\left({\frac {\mathbf {r} }{r}}\right).$

Для специального случая центральной силы, зависящей обратно пропорциональной квадрату расстояния $f(r)={\frac {-k}{r^{2}}}$, последнее выражение равно

${\frac {d}{dt}}(\mathbf {p} \times \mathbf {L} )=mk{\frac {d}{dt}}\left({\frac {\mathbf {r} }{r}}\right)={\frac {d}{dt}}(mk\mathbf {\hat {r}} ).$

Тогда $\mathbf {A}$ сохраняется в этом случае

${\frac {d}{dt}}\mathbf {A} ={\frac {d}{dt}}(\mathbf {p} \times \mathbf {L} )-{\frac {d}{dt}}(mk\mathbf {\hat {r}} )=0.$

Как показано ниже, вектор Лапласа — Рунге — Ленца $\mathbf {A}$ является частным случаем обобщённого сохраняющегося вектора ${\mathcal {A}}$, который может быть определён для любой центральной силы. Однако большинство центральных сил не формируют замкнутых орбит (см. теорема Бертрана), аналогичный вектор ${\mathcal {A}}$ редко имеет простое определение и в общем случае представляет собой многозначную функцию угла $\theta$ между $\mathbf {r}$ и ${\mathcal {A}}$.

Изменение под действием возмущающих центральных сил

Рис. 5. Медленно прецессирующая эллиптическая орбита, с эксцентриситетом $\scriptstyle e=0{,}9$. Такая прецессия возникает в проблеме Кеплера, если притягивающая центральная сила немного отличается от закона тяготения Ньютона. Скорость прецессии можно вычислить, используя приведённые в параграфе формулы.

Во многих практических проблемах, типа планетарного движения, взаимодействие между двумя телами только приблизительно зависит обратно пропорционально квадрату расстояния. В таких случаях вектор Лапласа — Рунге — Ленца $\mathbf {A}$ не постоянен. Однако, если возмущающий потенциал $h(r)$ зависит только от расстояния, то полная энергия $E$ и вектор углового момента $\mathbf {L}$ сохраняются. Поэтому траектория движения всё ещё находится в перпендикулярной к $\mathbf {L}$ плоскости, и величина $A$ сохраняется, согласно уравнению $A^{2}=m^{2}k^{2}+2mEL^{2}$. Следовательно, направление $\mathbf {A}$ медленно вращается по орбите в плоскости. Используя каноническую теорию возмущений и координаты действие-угол, можно прямо показать, что $\mathbf {A}$ вращается со скоростью

${\frac {\partial }{\partial L}}\langle h(r)\rangle ={\frac {\partial }{\partial L}}\left\{{\frac {1}{T}}\int \limits _{0}^{T}h(r)\,dt\right\}={\frac {\partial }{\partial L}}\left\{{\frac {m}{L^{2}}}\int \limits _{0}^{2\pi }r^{2}h(r)\,d\theta \right\},$

где $T$ — период орбитального движения и равенство $L\,dt=mr^{2}\,d\theta$ использовалось, чтобы преобразовать интеграл по времени в интеграл по углу (рис. 5). Например, принимая во внимание эффекты общей теории относительности, приходим к добавке, которая в отличие от обычной гравитационной силы Ньютона зависит обратно пропорционально кубу расстояния:

$h(r)={\frac {kL^{2}}{m^{2}c^{2}}}\left({\frac {1}{r^{3}}}\right).$

Подставляя эту функцию в интеграл и используя уравнение

${\frac {1}{r}}={\frac {mk}{L^{2}}}\left(1+{\frac {A}{mk}}\cos \theta \right),$

чтобы выразить $r$ в терминах $\theta$, скорость прецессии перицентра, вызванная этим возмущением, запишется в виде

${\frac {6\pi k^{2}}{TL^{2}c^{2}}}.$

которая близка по значению к величине прецессии для Меркурия необъяснённой ньютоновской теорией гравитации. Это выражение используется для оценки прецессии, связанной с поправками общей теории относительности для двойных пульсаров. Это согласие с экспериментом является сильным аргументом в пользу общей теории относительности.

Теория групп

Теорема Нётер

Теорема Нётер утверждает, что инфинитезимальная вариация обобщённых координат физической системы

$\delta q_{i}=\varepsilon g_{i}(\mathbf {q} ,\;\mathbf {\dot {q}} ,\;t)$

вызывает изменение функции Лагранжа в первом порядке на величину полной производной по времени

$\delta L=\varepsilon {\frac {d}{dt}}G(\mathbf {q} ,\;t)\,,$

что соответствует сохранению величины

$J=-G+\sum _{i}g_{i}\left({\frac {\partial L}{\partial {\dot {q}}_{i}}}\right)\,.$

Эта компонента вектора Лапласа — Рунге — Ленца $A_{s}$ соответствует вариации координат

$\delta _{s}x_{i}={\frac {\varepsilon }{2}}\left[2p_{i}x_{s}-x_{i}p_{s}-\delta _{is}\left(\mathbf {r} \cdot \mathbf {p} \right)\right],$

где $i$ равняется 1, 2 и 3, а $x_{i}$ и ${\dot {x}}_{i}$ — $i$-е компоненты векторов положения $\mathbf {r}$ и скорости $\mathbf {\dot {r}}$, соответственно. Функция Лагранжа данной системы

$L={\frac {m{\dot {r}}^{2}}{2}}+{\frac {k}{r}}.$

Получающееся изменение в первом порядке малости для функции Лагранжа запишется как

$\delta L={\frac {1}{2}}\varepsilon mk{\frac {d}{dt}}\left({\frac {x_{s}}{r}}\right).$

Это приводит к сохранению компоненты $A_{s}$

$A_{s}=\left[p^{2}x_{s}-p_{s}\ \left(\mathbf {r} \cdot \mathbf {p} \right)\right]-mk\left({\frac {x_{s}}{r}}\right)=\left[\mathbf {p} \times \left(\mathbf {r} \times \mathbf {p} \right)\right]_{s}-mk\left({\frac {x_{s}}{r}}\right).$

Преобразование Ли

Рис. 6. Преобразование Ли, из которого выводится сохранение вектора Лапласа — Рунге — Ленца $\scriptstyle \mathbf {A}$. Когда масштабируемый параметр $\scriptstyle \lambda$ изменяется, энергия и угловой момент тоже меняются, но эксцентриситет $\scriptstyle e$ и вектор $\scriptstyle \mathbf {A}$ не изменяются.

Существует другой метод вывода сохранения вектора Лапласа — Рунге — Ленца, использующий вариацию координат без привлечения скоростей. Масштабирование координат $\mathbf {r}$ и времени $t$ с разной степенью параметра $\lambda$ (рис. 6)

$t\to \lambda ^{3}t,\;\mathbf {r} \to \lambda ^{2}\mathbf {r} ,\;\mathbf {p} \to {\frac {1}{\lambda }}\mathbf {p}$

изменяет полный угловой момент $L$ и энергию $E$:

$L\to \lambda L,\;E\to {\frac {1}{\lambda ^{2}}}E$

— но сохраняет произведение $EL^{2}$. Отсюда следует, что эксцентриситет $e$ и величина $A$ сохраняются в уже упомянутом ранее уравнении

$A^{2}=m^{2}k^{2}e^{2}=m^{2}k^{2}+2mEL^{2}.$

Направление $\mathbf {A}$ также сохраняется, поскольку полуоси не изменяются при скалировании. Это преобразование оставляет верным третий закон Кеплера, а именно то, что полуось $a$ и период $T$ формируют константу $T^{2}/a^{3}$.

Скобки Пуассона

Для трёх компонент $L_{i}$ вектора углового момента $\mathbf {L}$ можно определить скобки Пуассона

$[L_{i},\;L_{j}]=-\sum _{s=1}^{3}\varepsilon _{ijs}L_{s},$

где индекс $i$ пробегает значения 1, 2, 3 и $\varepsilon _{ijs}$ — абсолютно антисимметричный тензор, то есть символ Леви-Чивита (третий индекс суммирования $s$, чтобы не путать с силовым параметром $k$, определённым выше). В качестве скобок Пуассона используются квадратные скобки (а не фигурные), как и в литературе и, в том числе, чтобы интерпретировать их как квантовомеханические коммутационные соотношения в следующем разделе.

Как показано выше, изменённый вектор Лапласа — Рунге — Ленца $\mathbf {D}$ можно определить с той же размерностью, что и угловой момент, разделив $\mathbf {A}$ на $p_{0}$. Скобка Пуассона $\mathbf {D}$ с вектором углового момента $\mathbf {L}$ запишется в похожем виде

$[D_{i},\;L_{j}]=\sum _{s=1}^{3}\varepsilon _{ijs}D_{s}.$

Скобка Пуассона $\mathbf {D}$ с $\mathbf {D}$ зависит от знака $E$, то есть когда полная энергия $E$ отрицательна (эллиптические орбиты под действием центральной силы, зависящей обратно пропорционально квадрату расстояния) или положительная (гиперболические орбиты). Для отрицательных энергий скобки Пуассона примут вид

$[D_{i},\;D_{j}]=-\sum _{s=1}^{3}\varepsilon _{ijs}L_{s}.$

В то время как для положительных энергий скобки Пуассона имеют противоположный знак

$[D_{i},\;D_{j}]=\sum _{s=1}^{3}\varepsilon _{ijs}L_{s}.$

Инварианты Казимира для отрицательных энергий определяются посредством следующих соотношений

$C_{1}=\mathbf {D} \cdot \mathbf {D} +\mathbf {L} \cdot \mathbf {L} ={\frac {mk^{2}}{2|E|}},$

$C_{2}=\mathbf {D} \cdot \mathbf {L} =0$

и мы имеем нулевые скобки Пуассона для всех компонент $\mathbf {D}$ и $\mathbf {L}$

$[C_{1},\;L_{i}]=[C_{1},\;D_{i}]=[C_{2},\;L_{i}]=[C_{2},\;D_{i}]=0.$

$C_{2}$ равен нулю, из-за ортогональности векторов. Однако другой инвариант $C_{1}$ нетривиален и зависит только от $m$, $k$ и $E$. Этот инвариант можно использовать для вывода спектра атома водорода, используя только квантовомеханическое каноническое коммутационное соотношение, вместо более сложного уравнения Шрёдингера.

Законы сохранения и симметрия

Вариация координаты приводит к сохранению длины вектора Лапласа — Рунге — Ленца (см. теорема Нётер). Это сохранение можно рассматривать как некоторую симметрию системы. В классической механике, симметрии — непрерывные операции, которые отображают одну орбиту на другую, не изменяя энергию системы; в квантовой механике, симметрии — непрерывные операции, которые смешивают атомные орбитали, не изменяя полную энергию. Например, любая центральная сила приводя к сохранению углового момента $\mathbf {L}$. В физике обычно встречаются консервативные центральные силы, обладающие симметрией группы вращения SO(3). Классически, полное вращение системы не затрагивает энергию орбиты; квантовомеханически, вращения смешивают сферические функции с тем же самым квантовым числом $l$ (вырожденные состояния), не изменяя энергию.

Рис. 7. Семейство кругов годографа импульса для заданной энергии $\scriptstyle l$. Все круги проходят через две точки $\scriptstyle \pm p_{0}=\pm {\sqrt {2m|E|}}$ на оси $\scriptstyle p_{x}$ (сравните с рис. 3). Это семейство годографов соответствует семейству окружностей Аполлония, и $\scriptstyle \sigma$ изоповерхностям биполярных координат.

Симметрия повышается для центральной силы, обратной квадрату расстояния. Специфическая симметрия проблемы Кеплера приводит к сохранению как вектора углового момента $\mathbf {L}$, так и вектора Лапласа — Рунге — Ленца $\mathbf {A}$ (как определено выше) и в квантовой механике гарантирует, что уровни энергии атома водорода не зависят от квантовых чисел углового момента $l$ и $m$. Симметрия является более тонкой, потому что операция симметрии должна осуществляться в пространстве большей размерности; такие симметрии часто называют скрытыми симметриями. С точки зрения классической механики более высокая симметрия задачи Кеплера учитывает непрерывные изменения орбит, которые сохраняют энергию, но не угловой момент. Другими словами, орбиты с одинаковой энергией, но различными угловыми моментами (эксцентриситетом) могут быть преобразованы непрерывно друг в друга. Квантовомеханически это соответствует смешиванию орбиталей, которые отличаются квантовыми числами $l$ и $m$, например, атомные орбитали типа $s$ ($l=0$) и $p$ ($l=1$). Такое смешивание нельзя произвести с обычными трёхмерными трансляциями или вращениями, но оно эквивалентно вращению в пространстве с более высоким измерением.

Связанная система с отрицательной полной энергией обладает симметрией SO(4), которая сохраняет длину четырёхмерных векторов

$|\mathbf {e} |^{2}=e_{1}^{2}+e_{2}^{2}+e_{3}^{2}+e_{4}^{2}\,.$

В 1935 году Владимир Фок показал, что квантовомеханическая проблема Кеплера эквивалентна проблеме свободной частицы, ограниченной четырёхмерной гиперсферой. В частности, Фок показал, что волновая функция уравнения Шрёдингера в пространстве импульсов для проблемы Кеплера представляет собой четырёхмерное обобщение стереографической проекции сферических функций из 3-сферы в трёхмерное пространство. Вращение гиперсферы и перепроектирование приводит к непрерывному преобразованию эллиптических орбит, не изменяющему энергию; квантовомеханически это соответствует смешиванию всех орбиталей с одинаковым главным квантовым числом $n$. Валентин Баргман отметил впоследствии, что скобки Пуассона для вектора углового момента $\mathbf {L}$ и скалированного вектора Лапласа — Рунге — Ленца $\mathbf {D}$ формируют алгебру Ли для $SO(4)$. Проще говоря, эти шесть величин $\mathbf {D}$ и $\mathbf {L}$ соответствуют шести сохраняющимся угловым импульсам в четырёх измерениях, связанных с шестью возможными простыми вращениями в этом пространстве (есть шесть способов выбрать две оси из четырёх). Это заключение не подразумевает, что наша Вселенная — четырёхмерная гиперсфера; это просто означает, что эта специфическая проблема физики (проблема двух тел для центральной силы, зависящей обратно квадрату расстояния) математически эквивалентна свободной частице на четырёхмерной гиперсфере.

Рассеянная система с положительной полной энергией обладает симметрией SO(3,1), которая сохраняет длину 4-вектора в пространстве с метрикой Минковского

$ds^{2}=e_{1}^{2}+e_{2}^{2}+e_{3}^{2}-e_{4}^{2}.$

Фок и Баргман рассмотрели как отрицательные, так и положительные энергии. Они также были рассмотрены энциклопедически Бендером и Ициксоном. Недавнее исследование Ефимова С. П. показало, что результат В. Фока переносится из искривлённого импульсного пространства в 4-х мерное координатное. При этом переход от четырёхмерных сферических функций в физическое трёхмерное пространство возникает просто при замене четвёртой «лишней» координаты на мнимый радиус вектор $\imath r$. Найденное координатное пространство оказывается в теории «ближе», чем искривлённое пространство В. Фока.

Симметрия вращений в четырёхмерном пространстве

Рис. 8. Годограф импульса на рис. 7 соответствует стереографической проекции больших кругов из четырёхмерной $\scriptstyle \eta$ сферы единичного радиуса. Все большие круги пересекают $\scriptstyle \eta _{x}$ ось, которая направлена перпендикулярно странице. Проекция из северного полюса (единичный вектор $\scriptstyle \mathbf {w}$) к ($\scriptstyle \eta _{x}$-$\scriptstyle \eta _{y}$) плоскости, как показано для пурпурного годографа пунктирной чёрной линией. Большой круг на широте $\scriptstyle \alpha$ соответствует эксцентриситету $\scriptstyle e=\sin \alpha$. Цвета больших кругов, показанных здесь, соответствуют цветам их годографов на рис. 7.

Связь между проблемой Кеплера и вращениями в четырёхмерном пространстве SO(4) можно достаточно просто визуализировать. Пусть в четырёхмерном пространстве заданы декартовы координаты, которые обозначены $(w,\;x,\;y,\;z)$, где $(x,\;y,\;z)$ представляют декартовы координаты обычного положения трёхмерного вектора $\mathbf {r}$. Трёхмерный вектор импульса $\mathbf {p}$ связан с четырёхмерным вектором ${\boldsymbol {\eta }}$ на четырёхмерной единичной сфере посредством

${\boldsymbol {\eta }}={\frac {p^{2}-p_{0}^{2}}{p^{2}+p_{0}^{2}}}\mathbf {\hat {w}} +{\frac {2p_{0}}{p^{2}+p_{0}^{2}}}\mathbf {p} ={\frac {mk-rpp_{0}}{mk}}\mathbf {\hat {w}} +{\frac {rp_{0}}{mk}}\mathbf {p} ,$

где $\mathbf {\hat {w}}$ — единичный вектор вдоль новой оси $w$. Поскольку ${\boldsymbol {\eta }}$ имеет только три независимые компоненты, то этот вектор можно обратить, получив выражение для $\mathbf {p}$. Например, для компоненты $x$

$p_{x}=p_{0}{\frac {\eta _{x}}{1-\eta _{w}}}$

и аналогично для $p_{y}$ и $p_{z}$. Другими словами, трёхмерный вектор $\mathbf {p}$ является стереографической проекцией четырёхмерного вектора ${\boldsymbol {\eta }}$, умноженному на $p_{0}$ (рис. 8).

Без потери общности, можно устранить нормальную вращательную симметрию, выбирая декартовы координаты, где ось $z$ направлена вдоль вектора углового момента $\mathbf {L}$, и годограф импульса расположен как показано на рис. 7, с центрами кругов на оси $y$. Так как движение происходит в плоскости, а $\mathbf {p}$ и $L$ ортогональны, $p_{z}=\eta _{z}=0$, и внимание можно сосредоточить на трёхмерном векторе ${\boldsymbol {\eta }}=(\eta _{w},\;\eta _{x},\;\eta _{y})$. Семейство окружностей Аполлония годографов импульса (рис. 7) соответствует множеству больших кругов на трёхмерной сфере ${\boldsymbol {\eta }}$, все из которых пересекают ось $\eta _{x}$ в этих двух фокусах $\eta _{x}=\pm 1$, соответствующих фокусам годографа импульса при $p_{x}=\pm p_{0}$. Большие круги связаны простым вращением вокруг оси $\eta _{x}$ (рис. 8). Эта вращательная симметрия преобразует все орбиты с той же самой энергией друг в друга; однако, такое вращение ортогонально к обычным трёхмерным вращениям, так как оно преобразует четвёртое измерение $\eta _{w}$. Эта более высокая симметрия характерна для проблемы Кеплера и соответствует сохранению вектора Лапласа — Рунге — Ленца.

Изящное решение для проблемы Кеплера с использованием переменных угол-действие можно получить, избавляясь от избыточной четырёхмерной координаты ${\boldsymbol {\eta }}$ и используя эллиптические цилиндрические координаты $(\alpha ,\;\beta ,\;\varphi )$

$\eta _{w}=\mathrm {cn} \,\alpha \,\mathrm {cn} \,\beta ,$

$\eta _{x}=\mathrm {sn} \,\alpha \,\mathrm {dn} \,\beta \cos \varphi ,$

$\eta _{y}=\mathrm {sn} \,\alpha \,\mathrm {dn} \,\beta \sin \varphi ,$

$\eta _{z}=\mathrm {dn} \,\alpha \,\mathrm {sn} \,\beta ,$

где используются эллиптические функции Якоби: $\mathrm {sn}$, $\mathrm {cn}$ и $\mathrm {dn}$.

Применение и обобщения

Квантовая механика атома водорода

Рис. 9. Уровни энергии водородного атома, предсказанные с использованием коммутационных соотношений углового момента и векторных операторов Лапласа — Рунге — Ленца; эти уровни энергии были проверены экспериментально.

Скобки Пуассона дают простой способ для квантования классической системы. Коммутатор двух квантовомеханических операторов равняется скобке Пуассона соответствующих классических переменных, умноженной на $i\hbar$. Выполняя это квантование и вычисляя собственные значения $C_{1}$ оператора Казимира для проблемы Кеплера, Вольфганг Паули вывел энергетический спектр водородоподобного атома (рис. 9) и, таким образом, его атомный эмиссионный спектр. Это изящное решение было получено до получения уравнения Шрёдингера.

Особенность квантовомеханического оператора для вектора Лапласа — Рунге — Ленца $\mathbf {A}$ заключается в том, что импульс и операторы углового момента не коммутируют друг с другом, следовательно, векторное произведение $\mathbf {p}$ и $\mathbf {L}$ должно быть определено тщательно. Как правило, операторы в декартовой системе координат $A_{s}$ определены с помощью симметризованного произведения

$A_{s}=-mk{\hat {r}}_{s}+{\frac {1}{2}}\sum _{i=1}^{3}\sum _{j=1}^{3}\varepsilon _{sij}(p_{i}l_{j}-l_{i}p_{j}),$

из которого определяются соответствующие лестничные операторы

$A_{0}=A_{3},$

$A_{\pm 1}=\mp {\frac {1}{\sqrt {2}}}(A_{1}\pm iA_{2}).$

Нормированный оператор первого инварианта Казимира может быть определён подобным образом

$C_{1}=-{\frac {mk^{2}}{2\hbar ^{2}}}H^{-1}-I,$

где $H^{-1}$ — оператор, обратный к оператору энергии (гамильтониан) и $I$ — единичный оператор. Применяя эти лестничные операторы к собственным состояниям $|lmn\rangle$ операторов полного углового момента, азимутального углового момента и энергии, можно показать, что собственные состояния первого оператора Казимира задаются формулой $n^{2}-1$. Следовательно, уровни энергии даются выражением

$E_{n}=-{\frac {mk^{2}}{2\hbar ^{2}n^{2}}},$

которое идентично формуле Ридберга для атома водорода (рис. 9).

Обобщение на другие потенциалы и СТО

Вектор Лапласа — Рунге — Ленца был обобщён на другие потенциалы и даже на специальную теорию относительности. Наиболее общую форму этого вектора можно записать в виде

${\mathcal {A}}=\left({\frac {\partial \xi }{\partial u}}\right)(\mathbf {p} \times \mathbf {L} )+\left[\xi -u\left({\frac {\partial \xi }{\partial u}}\right)\right]L^{2}\mathbf {\hat {r}} \,,$

где $u=1/r$ (см. теорема Бертрана) и $\xi =\cos \theta$, с углом $\theta$, определённым как

$\theta =L\int \limits ^{u}{\frac {du}{\sqrt {m^{2}c^{2}(\gamma ^{2}-1)-L^{2}u^{2}}}}\,.$

Здесь $\gamma$ — релятивистский фактор. Как и раньше, можно получить сохраняющийся вектор бинормали $\mathbf {B}$, взяв векторное произведение с сохраняющимся вектором углового момента

${\mathcal {B}}=\mathbf {L} \times {\mathcal {A}}\,.$

Эти два вектора можно соединить в сохраняющийся двухкомпонентный тензор $W$

${\mathcal {W}}=\alpha {\mathcal {A}}\otimes {\mathcal {A}}+\beta {\mathcal {B}}\otimes {\mathcal {B}}\,.$

Для примера вычислим вектор Лапласа — Рунге — Ленца для нерелятивистского изотропного гармонического осциллятора. Для центральной силы

$\mathbf {F} (r)=-k\mathbf {r} \,,$

вектор углового момента сохраняется, и поэтому движение происходит в плоскости. Сохраняющийся тензор можно записать в более простом виде:

$\mathbf {W} ={\frac {1}{2m}}\mathbf {p} \otimes \mathbf {p} +{\frac {k}{2}}\mathbf {r} \otimes \mathbf {r} \,,$

однако векторы $p$ и $r$ не перпендикулярны, как $A$ и $B$. Соответствующий вектор Лапласа — Рунге — Ленца имеет более сложную форму записи

$\mathbf {A} ={\frac {1}{\sqrt {mr^{2}\omega _{0}A-mr^{2}E+L^{2}}}}\{(\mathbf {p} \times \mathbf {L} )+(mr\omega _{0}A-mrE)\mathbf {\hat {r}} \}\,,$

где $\omega _{0}={\sqrt {\frac {k}{m}}}$ — частота осциллятора.

Ссылки

• Leach, P.G.L.; G.P. Flessas. Generalisations of the Laplace — Runge — Lenz vector (англ.) // J. Nonlinear Math. Phys. : journal. — 2003. — Vol. 10. — P. 340—423. Статья посвящена обобщению вектора Лапласа — Рунге — Ленца на потенциалы, отличные от кулоновского. arxiv.org Архивная копия от 12 августа 2020 на Wayback Machine